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1. INTRODUCTION 

The Digital image may be corrupted by different types of noise such as Gaussian, impulse, sparkle, 
and many others. Impulse noise usually corrupts the digital images due to errors produced from the noisy 
sensors or communication channel. It has a short duration and a sharp length. If impulse noise corrupts group 
of pixels in the original image, then they will be replaced by the impulse values, and the other pixels remains 
intact. As a result, the new image will have two types of pixels, the firt one is the noisy pixels n's that 
represent the impulse noise and the other type is the original pixels s's. Besides, there are two types of 
impulse noise, the first one is fixed valued impulse noise and the other one is random valued impulse noise. 
In the fixed type, the noise may take one or more value from the range of [0,255]. Example of this type is salt 
and pepper noise. In which, the salt noise takes a value equal to 255 and the pepper noise takes a value equal 
to 0. In random type, the noise may take any value in the dynamic range of [0,255]. In this paper, fixed 
valued impulse noise is investigated and the original image is artificially corrupted, either by one, two, or 
three fixed values. Impulse noise n(i, j ) located at i, j in the noisy image is modeled by a binary Bernoulli 
process b(ij) [1].In which b(i, j) equals to “1” with a probability p and equals to “0” with a probability 
1-p. Thus, the noisy pixel x(i, j) is expressed as: 


x(i, JJ=bG, j nG, j) + (1-bG , j))sG, j) (1) 


The challenge here is to restore the noisy image in such a way that preserve the fine details of the 
original image in the recovered version. The famous filter that used to suppress the impulse noise is the 
median filter [2]. Median filter replaces the tested pixel with the median value of a pre-defined filtering 
window centered at the tested pixel. The drawback of the median filter is that it replaces every pixel in the 
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image regardless if the tested pixel is noisy or not. Thus, many of the original pixels will be replaced by noisy 
values. As a result, the resultant restored image looks burred, and lose many original details. For enhancing 
the performance of the median filter, family of median-based filters are proposed. Among them are weighted 
median filter [3], and center weighted median filter [4]. In both, more weight is given to the pixels that are 
most likely original. Other median based filters are tri-state median filter [5], progressive-switching median 
filter [6], multi-state median filter [7], and pixel-wise MAD-based (PWMAD) filter [8]. 

These filters are decision-based filters used for preserving more image details by firstly detecting the 
noisy pixels. If the pixel is detected as noisy, then it is restored. To effectively restore the highly corrupted 
images, some novel switching-based filters have been proposed such as switching median filter with 
boundary discriminative noise detection BDND) [9], median-type noise detectors and detail-preserving 
regularization [10], and an iterative adaptive fuzzy filter based on two stages technique [11] and many others 
as [12]-[14].The authors in [15]-[17] proposed a new method based on the value of the noisy pixel and the 
pixels in the tested window.The drawback of these filters is that they provide rstored images with poor visual 
quality at high noise rates. Recently, many methods based on spare representation are used in literature as the 
ones mentioned in [18]-[21]. Unfortunately, these methods improve the results at high noise rates but at the 
expense of image details or computational complexity. In his paper, we assume that the fixed valued Impulse 
noise may take one, two, or three values each with different probability as: 


v1 withprobbiyit p1 


— v2 withprobbifit p2 
xí, f= : eg 
v3 withprobbifit p3 


s(i, j) withprobbifit - p (2) 


Where, pl+p2+p3=p. If the image is corrupted only by one value, then let p2=p3=0. If the image is corrupted 
by two values each with probability of 50%, then let p3=0. Thus, a new technique is used in this paper for 
detecting the noisy pixel based on the probability of existence of the tested pixel. More specifically, we 
should find the probability of existence of the tested pixel in order to be detected as original or noisy pixel. 
Any pixel is detected as noisy pixel if its probability of existence is more than a specific threshold. 
The proposed idea is unique and delivers superior and an outstanding results as shown in the following 
sections. The adopted idea does not depend on the surrounding pixels if they are original or not, and does not 
depend on the mean, median, or any other parameter for the image. It depends on independent factor which 
determines the presence of any value in the image. 


2. RESEARCH METHOD 
The proposed algorithm is divided into two parts the first part explains the detection process and the 
second part demonstrates the recovery process as shown below. 


2.1. Detection process 

From the literature, it is clear that many of the researchers make use the correlation between the 
pixels in the filtering window to detect the noisy pixels. In this paper, a new principle based on the 
probability of existence is proposed in the detection process. In the sense that, the repetitions of each pixel in 
the image is calculated, to help judging whether the tested pixel is an original pixel or not. To this end, we 
gather all the pixels that have equal values in separated groups. The number of the pixels in each group is the 
figure that we seek for. To illustrate the detection process, let X denotes to a noisy image of mxl size. The 
pixels that have density values equal to v; are collected in a separate group g;, i=1: K, where K is the number 
of the groups. The number of the pixels in the group g; is denoted by n;. Thus, the group g;, niis defined as: 


Gini=Wi,1Vi,2Vi,g -4 Vinil GB) 


Two matrices from (3) are defined; matrix V that contains one value from each group and matrix N that 
encloses the number n; of the corresponding value v;. Therefore, the matrices V and N are described as: 


V =[v, Vo V3 wee VK] (4) 


N=[M No Ng asese ng] (5) 
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The probability of existence p, of a value v,, 1 <r < K, is defined as: 


Hr (6) 


“Ixm 


P, 
Matrix P that includes the probability of existence for all image values is defined as: 

P=[P; P2 P3 «+e Px) (7) 
Define a new vector P to include all the probabilities that are greater than a threshold value py, as: 

PHD; » Pi>Pip:15iSM,M<K) (8) 
From (8) one can conclude that: 

P=1—Si0yPj = Dyas Pi =1-P a 


Where P represents the accumulative probability of all the probabilities that are less than py, and it is usually 
small value as shown in Table 1. Also, it is clear that, 


zM p; > MP, or! D; -MP >0 (10) 
1—P-—MP,, >0 (11) 
Pin <a (12) 


Let P =0 then the maximum bound for pm is less than z . Also, at highly corrupted images, the ratio P is 
M M 


M 
small because the noise rate or the probability of existence of the noisy pixel is high with respect to that of 
the original pixel, P <<, thus this ratio can be neglected and the value of pm becomes as: 


1 
0 < Pih < — 13 
th © m S 


From Table 1, one can observe that for a sample of 24 images, py, in any original image does not exceed 
than 4%. That means, we could not find a pixel in any original image with probability of existence equal to or 
more than 4%. It is clear from Figure 1 and Table 1, for Lena image, the pixel of the value 224 has the 
maximum probability of existence equal to 1.5%. Thus, the pixel in the noisy image that has probability of 
existence p, more than py, is considered noisy pixel, as described in Figure 2. Consequently, if p, > pm, then 
the pixel value v, at position r in matrix V is deemed noisy pixel as: 


n 
if = r 
1 B Tam. 


> Py, then vy is noisy vlue (14) 

xm 

At the end, the value v, is flagged as Z in the corresponding position in an image F of zero elements. As an 
example, noisy image corrupted with 80% salt and pepper noise rate has two noisy pixels each has 
probability of existence pr = 40% which is more than pm. Thus, we have only two values, i.e., M=2, that 
satisfy the condition in (14). During this paper, we consider py,= 5% which is less than 1/M, as indicated 
in 13 as: 


noise rate 
P — 


r 


> 5% <5 (15) 
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The following example represents the probability of existence for original and noisy values, let P and V as: 


P={ 0.03 0.03 0.02 0.02 0.4 0.03 0.03 0.02 0.02 0.4] 


V=[10 20 30 40 5 16 160 40 111 255] 


Then, it is clear that Ps=40%, Pj9=40%, and probability of original values is P=1-0.8=0.2 or 20%. 


number of each pixel in the image 


100 

50 

O Laie a hi hi 
50 100 


histogram of lena image 


u 


150 200 
pixel values in the image 


Figure 1. Illustrates the replication of each pixel in Lena image 


300 


Table 1. Maximum Probability of Existence in different known Images 


Image (value) ratio% Image (value) ratio% Image (value) ratio% 
Lena (224) 1.5 Lake (150) 1.31 Truck (129) 2.35 
Bridge (82) 0.87 Moon (134) 1.85 Tank (149) 2.91 
Pepper (199) 1.48 Airplane (207) 3.06 Bear (0) 3.05 
Boot (147) 2.0 Elaine (133) 0.90 Opera (144) 1.89 
Pentagon (130) 1.81 Gun (154) 1.42 Paper machine (206) 2.36 
Einstein (44 2.42 Water (106) 1.41 Couple (3) 3.38 
Watch (24) Sak Brandyrose (154) 2.28 Waterfall (117) 2.04 


Vie 


Image 


S VK 


“Jo ge(vkn) pr=nx/lm | 


q Bed | penn | 


v; is noisy pixel 


Figure 2. Illustrates the block diagram of the proposed detection technique 


2.2. Restoration process 


A pixel x, flagged as fj = 1 in F is replaced by the median value of its neighboring pixels Yj. Any 
pixel detected as noisy pixel and belong to Yj is excluded and the median is taken for the rest good pixels. 
This process is implemented recursively, in the sense that the current restored pixel is used in the restoration 
of the next noisy pixel in the subsequent windows. A small window W of 3x3 size is enough to be used in the 
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restoration process, but in the case of heavily corrupted image, noise rate > 70%, a window of 5x5 size or 
more is used. The neighboring pixels Y; and the restored version Xj rest Of the noisy pixel x; are defined as: 


Yi = { xi—s,j—t x fi—s,j—t | s,t eW}? 


$ 16 
Xij rest = median( Yj; ) te) 


The new method has the ability to detect the noisy pixel with zero error in large number of images. Thus, 
more powerful restoration process increases the performance of the proposed method and the visual quality 
of the recovered images. 


3. RESULTS AND ANALYSIS 

The proposed algorithm is compared with different known methods in terms of visual quality, 
processing time, and peak signal to noise ratio, PSNR. Note that the values of PSNR demonstrate the degree 
of sameness or closeness between the restored values and the original ones.Tested images are corrupted 
artificially with different noisy values. They are corrupted either with one, two, or three noisy values based 
on the MATLAB "imnoise" command. The strength of the proposed algorithm is showm clearly through its 
performance at high noise rates. Most of the known methods perform well at low noise rates, but the problem 
when the noise rate becomes > 40%. In this case, many of the well-known methods detect wrongly the true 
pixels. Therefore, maintaining the fine image details while detecting the noisy pixels is a challenge task. The 
advantage of the adopted algorithm is that it depends on a novel and unique feature explored during this 
research.This feature denotes that the structure and the view of any real image should have pixels its values 
fluctuate regularly and with some correlation.Based to this feature, the proposed method is capable of 
detecting all the noisy pixels correctly with a 100% accuracy, either at high or low noise rates and for large 
numbers of images. In the proposed approach, it is not necessary to classify the pixels as original, noisy, or 
most likely noisy through the detection and restoration processes as in many known methods. More 
specifically, the principle of the most known methods depends on one or more parameter or value such as the 
mean, median, standard deviation, noisy value, distance between the neighboring pixel, window size, 
correlation among the pixels, or to one or more threshold value. 

Unfortunately, these parameters or values fail to detect the noisy pixels, particurly at high noise 
rates. Because most of the pixels at this noise rates are noisy and therefore they involve in the detection 
process. Besides, it is difficult to estimate the original value of any pixel if most of their surrounding pixels 
are noisy. On the other hand, the principle of the new approach depends only on the weight of each pixel in 
the image.In other words, the assumed method indicates that the existance of any original pixel in any image 
is bounded and not exceeds a specific limit. But, the problem is how we can determine this limit to be a 
benchmark for detecting the noisy pixels. It is shown in Table land Figure 1 that the probability of finding 
any pixel in many original images is very small and not exceeds than 4%. Thus, the detection process states 
that if any pixel has probability of existence p, > threshold probability py, then it is considered noisy pixel. 
Threshold probability that used during the simulation is pm =5%. To demonstrate the power of the proposed 
method, extensive of experiments are performed and the results are shown numerically and visually through 
the following tables and figures. 

It is obvious from the results shown in Table 2, Table 3, and Table 4 that the proposed method 
achieves superior results comparing with the other known methods. It is also apparent from these tables that 
the new method provides the best PSNR values either at low or high noise rates in restoring different images 
corrupted with 0 and 255 noisy values. It is notable that, the performance of the proposed method in restoring 
corrupted Lena image is better than its performance in restoring bridge and baboon images. The reason is that 
the neighboring pixels in the bridge and baboon images are not similar enough to each other. As a result, the 
restored pixel during the recovery process does not reflect the optimum estimated value. Thus, the same 
detection method with more powerful recovery process would deliver better performance. Table 5 illustrates 
the performance of different methods in term of PSNR values when Lena image is corrupted by two different 
values 20 and 200. It is clear that the new method delivers superior results at low and high noise rates. 
Figure 3 compares the performance of the best two known methods with the proposed one in restoring 
corrupted Lena image with 90%, 70%, and 50% noise rates. Figure 4, demonstrations the performance of the 
proposed technique in restoring corrupted versions of boat, bridge, and baboon images. They corrupted at 
60% noise rate by the salt and pepper noise. It is notable that the proposed one outperforms clearly the others. 
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Table 2. Comparison for Restoring Corrupted Lena Image with 0 and 255 in Terms of PSNR (dB) 


Method/ 10% 20% 30% 40% 50% 60% 70% 80% 90% 

noise rate 
PSM [6] 3306 3000 2791 2620 2392 2163 1912 1548 1LI3 
ACWMF[4] 3721 34.75 3187 2854 2495 2062 1616 1232 825 
PWMAD[8] 3594 2659 2023 1605 1300 1064 873 718 589 
TSM [5] 25.55 2256 1867 1501 1214 982 795 651 543 
MSM [7] 38.24 3049 2475 1946 1550 1235 960 744 5.95 
NEW 38.34 35.65 3416 3300 32.013 3LI9 30.38 29.24 26.76 


Table 3. Comparison for Restoring Corrupted Bridge Image with 0 and 255 in Terms of PSNR (dB) 


Method/ 10% 20% 30% 40% 50% 60% 70% 80% 90% 

noise rate 
PSM [6] 28.86 26.75 24.76 2267 2088 1892 16.78 14.26 1045 
ACWMF [4] 27.27. 26.12 24.71 23.06 21.08 18.38 1548 1213 8.30 
PWMAD [8] 28.90 24.05 19.06 15.38 12.70 10.63 8.94 7.48 6.24 
TSM [5] 20.92 1995 17.50 14.76 12.18 9.99 8.23 6.85 5.84 
MSM [7] 28.36 26.20 23.01 1884 1543 12.26 9.70 7.10 6.31 
NEW 34.81 31.74 29.76 28.21 2681 25.45 23.97 22.25 19.97 


Table 4. Comparison for Restoring Corrupted Babbon Image with 0 and 255 in Terms of PSNR (dB) 


Method/ 10% 20% 30% 40% 50% 60% 70% 80% 90% 

noise rate 
PSM [6] 2417 2325 2246 2167 2072 1969 1823 +1621 1188 
ACWMF[4] 2373 2308 2232 2139 1986 1786 1514 1197 858 
PWMAD[8] 2569 2196 1791 1479 1231 1043 876 740 617 
TSM [5] 1972 1883 1674 1423 1183 981 816 682 58I 
MSM [7] 2452 2324 2110 1814 1493 1209 974 765 628 
NEW 3216 2910 27.15 25.74 2450 2340 2235 21.20 19.86 


Table 5. Comparison for Restoring Corrupted Lena Image with 20 and 200 in Terms of PSNR (dB) 


Method! 10% 20% 30% 40% 50% 60% 70% 80% 90% 

noise rate 
PSM [6] 33.16 3027 2818 2639 2410 2229 1989 1693 13.08 
ACWMF [4] 3652 3369 3039 2668 2306 1997 1681 1325 979 
PWMAD[8] 3666 2896 2238 1802 1492 1228 1026 865 727 
TSM [5] 2555 2339 1970 1619 1344 1122 946 807 698 
MSM [7] 36.53 3148 2540 2050 1695 1351 11.00 894 7234 
NEW 4277 39.47 3727 3539 3355 3179 3054 2932 2690 


Figure 3. Comparison between existing methods and the proposed one for restoring salt and pepper noise 
[0,255] corrupted Lena images: (a) New (b) ACWM (c) PSM (d) corrupted Lena images. 1“, 2" and 3" 


images in row (d) are corrupted by 90%, 70, and 50% noise rates, respectively 
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Figure 4. Comparison between existing methods and the proposed one for restoring 60% salt and pepper 
noise [0,255] corrupted three differnt images. Corrupted images: (a) New (b) ACWM (c) PSM 
(d) corrupted boat, bridge, and baboon images 


Figure 5 illustrates the performance of the proposed method in the restoration of 60% corrupted 
pepper images either by one, two or three noisy values. It is clear that the proposed method shows superior 
performance than the others, particularly at 70% and 90% noise rates. It is also clear that there are residual 
noise and blur in the restored images taken from the other methods. Moreover, one can note that the 
increment in the noisy pixels does not affect the performance of the proposed method. The only thing that 
shakes its performance is the efficiency of the restoration process in replacing the noisy pixel by an estimated 
value similar to the original one. Thus, in the restoration process, it is recommended to use a wider filtering 
window at high noise rates to make sure that the filtering window has an enough number of original or 
already restored pixels for well-estimate the noisy pixel. Table 6 shows the performance of the new filter in 
restoring Lena image corrupted with one noisy value (30) at different noisy rates. Table 7 illustrates the 
performance of different methods in restoring lina image corrupted with three noisy values 0, 255, and 30 
with rates of 50%, 25% and 25%, respectively. It is clear that the proposed method is an extraordinary 
approach in terms of PSNR, and particularly at high noise rates. 

The other known methods have difficulty to detect the noisy pixels, mostly at highly corrupted 
images. The reason is that the other methods almost depend on threshold values based on different 
parameters such as mean, median, or standard deviation of the surrounding pixels. These parameters almost 
provide inaccurate values when most of the surrounding pixels are noisy. Figure 6 depicts the consumed time 
taken by the PSM, ACWMF, and the proposed one in restoring pepper image at different noise rates. It is 
clear that, the new method outperforms the others in terms of the processing time in seconds.Thus, 
the proposed filter is simple, converges very fast, and delivers superior results. MATLAB program, CPU 
of 1.73 GHz, and RAM of 1GB are used in all the simulated experiments. 
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Figure 5. Comparison between existing methods and the proposed one for restoring corrupted versions of 
pepper image: row (a) New, row (b) ACWM, row (c) PSM, row (d) corrupted versions of pepper images by 
60% noise rates. In row (d) first pepper image corrupted by one value [30], 2" image corrupted by two 

values [20,200], and 3" image corrupted by three values [0,255,30], respectively 


Table 6. Comparison for Corrupted Lena Image by one value “30” in Terms of PSNR (dB) 


Method/ 


Hoete 10% 20% 30% 40% 50% 60% 70% 80% 90% 
PSM [6] 37.53 33.93 29.43 20.09 6.97 4.19 4.15 4.15 4.15 
ACWMF [4] 36.11 30.96 23.15 14.95 7.28 4.51 4.19 4.15 4.15 
PWMAD [8] 31.90 20.09 13.63 9.63 7.11 5.52 4.64 4.24 4.15 
TSM [5] 23.77 17.29 12.38 9.21 7.07 5.67 4.79 4.32 4.16 
MSM [7] 33.50 21.95 14.88 10.14 7.05 5.32 4.50 4.21 4.15 
NEW 42.99 39.56 37.28 35.40 33.66 31.88 30.65 29.38 26.78 


Table 7. Comparison for Corrupted Lena Image by 0 (50%),255 (25%), and 30 (25%), in Terms 


of PSNR (dB) 
Method/ 10% 20% 30% 40% 50% 60% 70% 80% 90% 
noise rate 
PSM [6] 34.68 31.93 29.66 27.16 22.65 14.37 4.56 3.74 3.69 
ACWMF[4] 36.93 33.91 28.96 23.25 17.05 10.68 5.81 4.46 4.18 
PWMAD [8] 35.23 24.38 17.59 13.12 10.11 7.96 6.35 5.24 4.49 
TSM [5] 25.29 21.41 16.62 12.85 10.12 8.03 6.49 5.37 4.65 
MSM [7] 36.43 28.71 21.78 16.43 12.14 8.82 6.40 5.08 4.48 
NEW 39.06 37.00 35.21 3343 3166 2961 3046 29.28 26.67 
COMPARISON IN TERMS OF COMPUTIONAL COMPLEXITY 
300 T T T T T 
3 250 A 
2 = ase 
E3 150 
È i i a. = 
2 o l 
r e a 7 e 
10 20 30 40 60 70 80 


Figure 6. Comparison between different known methods and the new one in terms of computional 
complexity due to restoring pepper image corrupted at different noise rates 


impulse noise rate 
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4. CONCLUSION 

In this paper, a new filter based on a unique idea is proposed. It has been shown that the presence of 
any pixel in many tested images is bounded. Thus, the tested pixel is considered noisy pixel if its probability 
of existence is more than specific threshold. Extensive simulation results prove that the proposed method 
provide outstanding results at low and high noise rates. The reason is that its detection process does not 
depend on any correlation between the neighboring pixels which vary from window to window and from 
region to region. Moreover, the new filter is very fast and easy to implement. 
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